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Abstract 

It is now clear that heterogeneous reactions play an impor- 
tant role in controlling the concentration of ozone in the 
stratosphere. In this work, the loss of N 2 O 5 on ice sub- 
strates has been studied in a flow reactor in an attempt to 
gain a more fundamental understanding of these reactions. 
The apparent reaction probability in this system was found 
to decrease as the substrate was exposed to N 2 O 5 . A model 
which corrected for the loss of surface sites was developed 
and although it appears to fit the data for a given exper- 
iment quite well, it is concluded that the loss of reactive 
sites is not the full explanation. In addition, the results of 
an experimental and modelling study suggest that reaction 
on the internal surface of the ice substrates is not a major 
loss mechanism for N 2 Os in the current work. 

1 Introduction 

It is now recognized that reactions which occur on the 
surfaces of aerosols in the stratosphere have a significant 
effect on the concentrations of ozone. Before laboratory 
measurements can be used with confidence in atmospheric 
models to explain and predict the effects of heterogeneous 
chemistry in the stratosphere, much more detailed informa- 
tion about these processes is required. In this study, the 
reaction of N 2 O 5 with water ice substrates has been chosen 
as a sample system in order to address several fundamental 
questions about heterogeneous processes as studied in the 
laboratory. 

This system has been studied by several groups [Hanson 
and Ravishankara , 1991; Leu, 1988; Quinlan et a/., 1990; 
Tolbert, Rossi and Golden , 1988] and it was shown that the 
product of the reaction is nitric acid which is retained in 
the solid-phase. 

N 2 0 5( ,) + H 2 0 (#) — > 2HN0 3W (1) 


acid which is thought to react further with water to form 
nitric acid trihydrate (NAT) [Hanson and Ravishankara , 
1991; Kenner, Ryan and Plumb , 1992]. Using low initial 
concentrations of N 2 O 5 and short measurement times to 
minimize the build-up of HNO 3 product on the surface 
of the ice, Hanson and Ravishankara [Hanson and Rav- 
ishankara , 1991] estimated a value of «2xl0 -2 for the col- 
lisional reaction probability of N 2 O 5 on ice at 200 K. They 
found that the apparent reaction probability decreased on 
continued exposure of the ice to N 2 O 5 to a limiting value 
of 6xl0 -4 . They interpreted the latter value as being the 
reaction probability on NAT which they reasoned would 
have covered the surface completely at the doses of N 2 O 5 
used. 

If the observed decrease in reaction probability is only 
due to the gradual conversion of the original ice surface 
into one which is covered with NAT, it should be possible 
to model the experiments and obtain a value of the reaction 
probability of N 2 0 5 on pure ice, 7 i ce . Then using the model 
it should be possible to predict the reaction probability of 
N 2 O 5 on surfaces which are intermediate between pure ice 
and pure NAT. In this work the reaction probability of 
N 2 O 5 on ice substrates with submonolayer coverage of the 
surface by HNO3 was measured and modelled. 

The second question addressed in this work is the de- 
termination of the correct surface area to use in converting 
the observed wall-loss rate to a reaction probability. Most 
authors have assumed that the surface of their substrates 
are smooth and non-porous and have used a surface area 
calculated from the dimensions of their apparatus. Re- 
cently, Leu and coworkers [Keyser, Moore and Leu , 1991; 
Leu, Moore and Keyser , 1991] have observed that, for ni- 
tric acid hydrate substrates, the measured wall-loss rate 
increases with the thickness of the film in a manner which 
suggests that reaction on internal surfaces of those films is 
important. They showed that correction factors as large as 
50 may need to be applied to laboratory data to account 
for the increased area of the substrate due to its porous 
structure. 

In this work, two approaches have been taken to de- 
termine the possible importance of reaction on internal 


where (g) and (s) refer to species in the gas and solid 
phases respectively. As the reaction proceeds, the surface 
of the substrate is increasingly contaminated with nitric 
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surfaces of the substrate in the current system. (1) The 
reaction probability of N 2 0s on ice substrates has been 
measured as a function of the thickness of the substrate. (2) 
A model has been developed to fit the data which includes 
diffusion into and reaction on the surface of idealized pores. 


2 Results and Discussion 

2.1 Experiments and Models Concern- 
ing the Conversion of Reactive Sites 

The flow reactor and detection system used in these 
studies have been described previously [Plumb and Ryan , 
1981; Ryan and Plumb , 1982] and its modification and ap- 
plication to the study of heterogeneous reactions at low 
temperatures is described elsewhere [Kenner , Ryan and 
Plumb , 1992], N 2 0s was prepared by a method similar to 
that described by Davidson et ai [Davidson et al. , 1978] 
and Cantrell et al\Cantrell et al. , 1987]and was detected 
at its NOj (m/e 46) fragment. 

Results for two experiments measuring the loss of N 2 0s 
on ice-coated walls of the flow reactor are shown in Figure 1 
as points in semilogarithmic plots of [N 2 Os] as a function of 
inlet distance. The error bars are the statistical uncertain- 
ties in each point. The points in Figure 1 were obtained by 
measuring [N 2 Os] at fixed positions 100 mm apart as the 
movable inlet is moved from an inlet distance of ~50 mm 
to an inlet distance of 450 mm and back i.e. the points 
were measured in order going round the loop in a counter- 
clockwise direction. The curves shown in Figure 1 are the 
results of fits using a model to be discussed below. 

In all experiments in the current investigation, loops 
such as those shown in Figure 1 were obtained and show 
that the reaction probability of N 2 0 5 on the surface changes 
during the course of an experiment. The shapes of the loops 
depend on a number of experimental variables. For exam- 
ple, with all other parameters held constant, the loops are 
broader with a more horizontal major axis for: i) higher 
[N 2 0 5 ]o (compare curves a and b in Figure 1) or ii) longer 
count times. 

For the purposes of modelling, it is possible to represent 
the overall process by a single reaction 

N 2 0 5(3) + s/H 2 0 (j) — > 2(HN0 3 ■ «II 2 0) h ( 2) 

For HN0 3 -nH 2 0( s ) — NAT, n — .‘3 and s/, the stoichiomet- 
ric factor, would be 7 if all water molecules of hydration 
are derived from the surface monolayer. A model has been 
constructed to fit the current experiments based on Reac- 
tion (2) in which sf and 7 are adjustable parameters whose 
values were obtained by fitting the experimental points us- 
ing non-linear least-squares fitting techniques based on the 
Levenberg-Marquardt method [Press ct al. , 1988]. This 
method includes uncertainty estimates for the fitted pa- 
rameters. 

As the experiment proceeds, the ice sites are lost due to 
the formation of nitric acid hydrates and the model must 



Inlet Distance (mm) 

Figure 1: N 2 0& concentration as a function of inlet dis- 
tance for ice-coated flow tube walls at a temperature of 
175 K for two different starting concentrations of N 2 05. 
The points are experimental data (converted to absolute 
concentrations using the experimentally determined sensi- 
tivity for N 2 Os) and the curves are the best fit results of 
the basic model described in the text. 


keep track of the surface concentration of ice and NAT sites 
at each position along the flow' tube as a function of time. 
The reaction probability at each point is assumed to be 
proportional to the number of ice sites, [H 2 0( s )(z,t)], 


7(z,t) 


[H 2 0 (s) (z,t)] 

[H 2 0 (s) ]o 


(3) 


where [H 2 0( a )]o is the surface concentration of ice sites on 
pure ice. The rate equations are integrated to determine 
the concentrations of N 2 Os and H 2 0( 5 ) as functions of posi- 
tion and time. The average concentration of N 2 ()5 surviv- 
ing to reach the detector is calculated for each position of 
the movable inlet . By setting the movable inlet positions, 
dwell time, flow velocity, inlet scan speed and [N ^Or. ]u in 
the model to match those of a particular experiment, it is 
possible to simulate that experiment. 

The surface concentration of ice sites at the beginning 
of the experiment was assumed to be IxiO 1 ’ cm' 2 . Be- 
cause of the way tin* model is constructed, the actual value 
chosen for the siirlace concentration has no effect, on the 
value of 7 obtained. However, the value ol *J is directly 
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Ice Thickness (/xm) 

Figure 2: The values of 7 from the basic model as a function of the thickness of the ice film. 


proportional to the choice of the initial surface concentra- 
tion of reactive sites. Based on the results of consecutive 
experiments performed using the same substrate, regener- 
ation of the ice sites through solid state diffusion and/or 
evaporation of the HNO 3 is considered to be sufficiently 
slow that it can be ignored for T<200 K. 

The basic model fits the data quite well as can be seen 
in Figure 1 , where the curves are the best fits to the data. 
From the good fits to the data obtained using the basic 
model, one can conclude that the main reason for the de- 
crease in the reactivity during the experiments is the grad- 
ual conversion of the more reactive ice sites into less reac- 
tive NAT sites. However, one also finds that experiments 
that differ only in the amount of exposure of the substrate 
to N 2 0& yield different values for the reaction probability 
on pure ice. This is a clear indication that the conversion 
of sites is not the full explanation. 

Several additional mechanisms have been considered, 
for example, higher than first-order dependence on [H^O^)], 
sites of different reactivities and involvement of more than 
just the top monolayer of the substrate. Models to fit the 
data have been developed based on each of these and al- 
though they each improve the fits to the data and reduce 
the dependence of 7 on t he experimental conditions, they 
introduce too many parameters about which no informa- 
tion is available and no conclusion can be reached as to the 
correct ness of any given model. 


2.2 Dependence of the Reaction Proba- 
bility on the Thickness of Ice and the 
Porous Substrate Model 

The dependence of the reaction probability of N 2 Os 
on ice has been investigated in this work for T = 190 K, 
bath gas pressure = 0.23 Torr and [N20s]o = 3.5xl0 10 
cm -3 for ice thicknesses, h,- ce > in the range 0<h tce <124 /im. 
The results are shown in Figure 2 where 7 mo d given by the 
basic model are plotted as a function of h lce . It will be 
seen that 7 mo d increases with h,- ce for 0<h, ce <7 fim and is 
independent of h tce for h tce >7 ^m. The authors believe 
that the lower values of 7 observed for ice thicknesses less 
than 7 are caused by the walls of the flow tube being 
incompletely covered with ice in those experiments. This 
view is supported by the large values of sf obtained for 
the experiments on thin ice substrates [A'enner, Ryan and 
Plumb , 1992]. 

The possible porosity of the ice in the current experi- 
ments has been tested further by developing a model which 
includes idealized pores to investigate what ranges of pore 
diameters and porosities (the ratio of the volume of internal 
space to the volume of solid) are capable of explaining the 
current data. By comparing the calculated loops with the 
experimental data for a range of values of the mean pore 
diameter and of the porosity of the substrate, it has been 
found that mean pore diameters of 0.1 /mi and porosities 
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of 0.03 are required to reproduce the data. These values 
of mean pore diameter and porosity were tested by fitting 
the same data sets used in evaluating the other models. 
The deduced values of 7 differed by <20% from those ob- 
tained using the basic model. The average value of \ 2 and 
the dependence of the calculated reaction probability on 
experimental conditions were similar to those obtained us- 
ing the other ‘improved’ models. Based on these results, 
it is concluded that in the current system reaction on the 
internal surface of the substrate is not a large effect but 
may explain some details of the current experiments. Any 
further evaluation of this model must await experimental 
measurements of the pore diameters and porosities of the 
substrates. 
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